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Repeated closed-loop control operations acting as piecewise-constant Liouville superoperators 
conditioned on the outcomes of regularly performed measurements may effectively be described by 
a fixed-point iteration for the density matrix. Even when all Liouville superoperators point to the 
completely mixed state, feedback of the measurement result may lead to a pure state, which can be 
interpreted as selective dampening of undesired states. Using a microscopic model, we exemplify this 
for a single qubit, which can be purified in an arbitrary single-qubit state by tuning the measurement 
direction and two qubits that may be purified towards a Bell state by applying a special continuous 
two- local measurement. The method does not require precise knowledge of decoherence channels 
and works for large reservoir temperatures provided measurement, processing, and control can be 
implemented in a continuous fashion. 
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PACS numbers: 03.65.Ta 03.65.Yz 03.67.Bg 03.67.Pp 

In quantum systems, one typically aims at avoiding 
decoherence that is often seen as arch-enemy of quan- 
tum computation [l|. Simply performing the computa- 
tion fast enough will fail for most applications as a size- 
scalable quantum system with large coherence times [2| 
is yet to be found. Therefore, advanced schemes have 
been proposed to reduce or inhibit decoherence in quan- 
tum computers such as e.g. quantum error correction ^j 
the use of decoherence-free subspaces Q, or open- loop 
control pH8!|. These schemes require quite sophisticated 
techniques which themselves may not be entirely robust 
against decoherence and/or control errors. 

The idea to revert the perspective by using decoherence 
constructively is not new |9- 12] but has recently again 
gained a lot of attention |13l - fl7| : The simplest exam- 
ple is to relax into the (possibly entangled) ground state 
of a defined system Hamiltonian, which only requires 
weak unspecific couplings to a low-temperature reser- 
voir with a sufficiently large energy gap above the ground 
AE ^ ksT. The limitation to the ground state can be 
overcome by using multiple reservoirs at different ther- 
mal equilibria and different chemical potentials. Then 
however, the general dissipative engineering paradigm re- 
quires to design the interactions such that the system is 
driven to the desired target state. For a general pure 
state to be stabilized, this problem is hard to solve, the 
solution is specific to the target state and may prove dif- 
ficult to implement experimentally. 

Since the success of the centrifugal governor used by 
James Watt to regulate the speed of steam engines, feed- 
back (closed-loop) control has found wide-spread applica- 
tion with todays standard examples including e.g. ther- 
mostats and automatic speed control in cars. Nowadays, 
the built-in mechanical self-regulation of the centrifugal 
governor is often replaced by electronic signal process- 
ing, which enables one to change the feedback protocol 
easily. The implementation of quantum feedback control 



requires the inclusion of the quantum-mechanical mea- 
surement process and to make the quantum-mechanical 
evolution after each measurement subject to control. In 
the conventional scheme [la, [l9|, quantum jumps de- 
tected during the time-interval At are fed back as an 
instantaneous (5-pulse into the system Hamiltonian. In 
the continuum measurement limit At — >■ 0, an effective 
master equation emerges, where the control operation ap- 
pears as an instantaneous unitary rotation of the density 
matrix [13-113 ■ 

In contrast, here we consider a general scheme of pe- 
riodic measurements at intervals At and discuss in de- 
tail the opposite case where the collapse due to projec- 
tive measurements is assumed to be much faster than 
At throughout. The scheme allows for weak measure- 
ments and in principle arbitrary control pulse sequences, 
but we will for simplicity of interpretation specialize on 
a piecewise-constant time-dependence of the control pa- 
rameters. In this scenario, even for Ai — >■ a master 
equation description does not typically emerge. The pur- 
pose of this paper is to show that, provided with dissipa- 
tion that drags the system only to the completely mixed 
state but allows for strength-tunable system-reservoir in- 
teractions, feedback of information obtained by measure- 
ments can be used to achieve purification and thereby 
puts even "clueless" dissipation to a productive role. 



PROPAGATION UNDER FEEDBACK 



A general quantum measurement is described by a set 
of measurement operators {Mm} satisfying the complete- 
ness relation J2m ^m^-m = 1 [H- Orthogonal projection 
operators MnMm = SnmMn are a well-known standard 
case. Under measurement outcome to, the density ma- 
trix becomes (we assume an instantaneous collapse of the 
wave function) 
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MmpMl 

P{m) 



(1) 



and the outcome probability is given by P{m) — 
Tr {M^Mmp}. Arranging the N^ elements of the den- 
sity matrix in a vector, this becomes 



P{m) 



MmP, 



(2) 



where A^,„ is a superoperator having an N'^ x N'^ matrix 
representation. In the following, we will use calligraphic 
symbols to denote superoperators. In principle, the mea- 
surement operators may depend also on the time interval 
At between two measurements, see appendix [Xj 

The most general evolution of the density matrix (in- 
cluding measurements, unitary and dissipative evolution) 
is governed by a trace- and positivity preserving map 

Pit + At) = Y,BaiAt)p{t)Bi{At) 

a 

= J2B^{At)p{t)^B{At)pit) (3) 

a 

with Kraus operators obeying Y^^Bl^{At)Ba{At) = 1 
and the corresponding superoperator B{At). Formally, 
feedback can be conveniently described by making the 
general evolution between two measurements dependent 
on the outcome of the previous measurement. For ex- 
ample, given density matrix p{t) right before one mea- 
surement, the density matrix conditioned on outcome m 
at the first measurement becomes right before the next 
measurement 



M^p{t)Ml ^^ 

Bm{At)MraP{t) 



p{t + At) n Y.B^.„,{Aty-^^^f^^Bl,^{At) 



P{m) 



(4) 



For an arbitrarily chosen observable, we perform a 
weighted average for its expectation value at time t + At 
over all measurement outcomes and conditioned evolu- 
tions 



(At+At) = Ti- { A 



J2>3,n{At)Mr 



Pit) , (5) 



where we have used superoperators for the ease of no- 
tation with Ap=Ap and the trace is generally mapped 
to multiplication with a row vector containing N en- 
tries of value 1 at positions containing populations and 
7V^ — N entries of value elsewhere. Demanding 

(A-t+At) = Tr {Ap{t + At)} defines an effective propaga- 
tor for repeated measurements and feedback 



VM^t)=Y,Bm{At)Mr 



(6) 



The above propagator may prove useful once microscopic 
parameters are linked to the superoperators Bm{At), 
which is however in prin ciple possible also for non- 
Markovian systems |23l - l25| . In appendix]^ we show that 



it yields an effective master equation as in the conven- 
tional Wiseman-Milburn scheme for non-projective quan- 
tum jump detection, unitary control and in the contin- 
uum limit At — > 0. 

In the following however, we will consider a Lindblad 
evolution with a single constant Lindblad superoperator 
Cm (wh ich may however also include Hamiltonian con- 
trol [2a|) between the measurements Bm{At) = e 
and time-independent measurement superoperators 



C^At 



VesiAt) 



^e- 



At 



Mr 



(7) 



Note that in contrast to the Wiseman-Milburn scheme 
(see also appendix IX)) . the time- independence of the mea- 
surement superoperators requires the measurement to be 
much faster than At throughout. The conditioning of the 
Cm on the previous measurement result m defines the 
feedback protocol. Microscopically, this dependence on 
the measurement result can be implemented by triggering 
switches of parameters in the Hamiltonian. We explicitly 
include the possibility that the stationary state (defined 
by CmP = 0) of each Liouvillian may be the completely 
mixed state with Tr{p^} = 1/A^. In the following, we 
will discuss the implications of this effective propagator 
for the control of qubits. 

Typically, decoherence in an open quantum system is 
modeled by the decay of off-diagonal matrix elements of 
the density matrix (coherences). In an appropriate ba- 
sis, this is formally represented by a block form of the 
Liouvillians £m, which leads to a decoupling of coher- 
ences and diagonal entries (populations). A prototypical 
example for this is the Born-Markov-secular approxima- 
tion master equation in the case of a non-degenerate sys- 
tem Hamiltonian [23| , where the block structure becomes 
explicit in the system energy eigenbasis. However, for 
quantum measurements the measurement superoperators 
Mm will not generally have the same block structure as 
the Liouvillian. It may therefore be expected that the 
effective propagator in Eq. ([7]) couples coherences and 
populations and thereby may lead to nontrivial effects. 

In spite of the completeness relation in the original 
Hilbert space, the measurement superoperators do not 
sum up to the identity X)m-^™ t^ ■'■• This implies 
that even in the continuum limit, where measurement, 
processing, and control are repeated at infinitesimally 
small time intervals At — ^ 0, the iteration equation 
p{t + At) = VcsiAt) p{t) with ([7]) does not converge to 
an effective master equation p ^ Cc^p- A prominent 
example for this is the quantum Zeno effect: For a closed 
system without measurements and without feedback, 
the action of the Liouvillian without measurement 
simply becomes p{t + At) = e^°^^p{t)=e-'"^*pe+'"^\ 
which just encodes the usual unitary evolution. 
With measurements and without feedback, the ef- 
fective propagator in Eq. ([7]) leads to p{t + At) = 



CoAt 



Em Mm] p(i)=e-'^^* [Em MmP{t)Ml] 



,+iHAt 



where the measurement (super-) operators may account 
for the freezing of the observed quantum state when 



At — >■ 0, as will be discussed in the next section in 
greater detail. 

A continuous description by means of a differential 
equation for expectation values may however still be pos- 
sible when the expectation values are consistent with the 
measurement operators, as will be shown below. 



A. Thermal reservoirs 

Let us first consider the simplest case where neither 
measurement nor control are applied: The continuous 
action of either Liouvillian would lead to an exponential 
decay of coherences |poi| {t) = e~^±*|poi| (0) with the 
dephasing rate 



II. DISSIPATIVE PURIFICATION OF A QUBIT 

As a first example we consider a single qubit pa- 
rameterized by Pauli matrices. At time intervals At, 
we perform strong projective measurements of a'^ de- 
scribed by the measurement operators M± = ^ [1 ± a^] 
with superoperator equivalents acting on the vector 
(poo , Pii 1 Poi , Pw) having the representation 



A|{4cos2(0±)7o 



-.)b+n+l-n]} ■ (10) 



M- 



1 1 ±1 ±1 

1 / 1 1 ±1 ±1 

4 I ±1 ±1 1 1 

±1 ±1 1 1 



(8) 



The feedback protocol is defined by applying an outcome- 
dependent Liouvillian derived microscopically from the 
system (S), bath (B) and interaction (SB) Hamiltonians 



-ff, 



(±) _ 



n 



^a^ 4±)=^c.46., 



Hg'> = A± (n± -CTJ^Y. ('^fe^fe + ''l^i) ' (9) 



the 



unit vectors n± — 
^±)) character- 



where A± > and 
(sin(6'±) cos{(j)±), sin(6'±) sin((/)±), cos( 
ize strength and direction (e.g., purely dephasing for 
n = Ez) of the dissipation, respectively. The bk are 
bosonic annihilation operators of a bath assumed at 
thermal equilibrium throughout. The case X± = X and 
n± = n corresponds to repeated measurements without 
feedback applied to an open quantum system. Applying 
the Born, Markov and secular approximations, the 
corresponding conditioned Liouville superoperators C± 
do not depend on (/)± and both have block structure in 
the system energy eigenbasis a^ |0/1) = (—1)'^°/-'^) |0/1), 
leading to a decoupling of coherences and populations, 
see appendix IBJ 

The ability to change the system-reservoir coupling will 
of course depend on the physical implementation of the 
qubit. In electronic setups for example (charge qubits), 
nearby quantum point contacts may not only function 
as detector: Changing the bias voltage for circuits in the 
vicinity of the charge qubit will induce a changed system- 
reservoir coupling. If in contrast the qubit is realized as 
the lowest two modes of a cavity, tuning the permeability 
of the cavity walls may yield a similar effect. 



Here, 7cu = / {e+'"^'' Be-^'"^^ B) e+^'^^dr is the Fourier 
transform of the bath correlation function with the 



Ek{f 



hhbk + h*i.b,. 1 when 



k^k 



coupling operator B 

thermal state at inverse temperature (3 is assumed. 
Analytically continuing the spectral coupling density 
J(a;) = 27r^^ \hk\ 5{u! — uik) to negative lo via 
J(— w) = —J{+co) one obtains with the Bose distri- 
bution 71b ('^) = [e^^"^ — l] the simple expression 
7lj — J(y>) [1 + n-B,{u!)]. For both Liouvillians, the pop- 
ulations will approach the thermalizcd stationary state 
characterized by a^ = (7_o - 7+n)/(7-o + 7-i-n)- Due 
to micro-reversibility, the correlation function obeys the 
Kubo-Martin-Schwinger [23] condition 7_n — e^'^^7+n. 

Therefore, for low temperatures fceT <C il, the qubit 
will simply decay towards its ground state, which is 
a trivial example of dissipation-induced purification in 
the eigenbasis of the system Hamiltonian. Note how- 
ever that the validity of the Markovian approximation is 
only expected when the coupling strength is smaller than 
the temperature, such that the relaxation speed may be 
slowed by scaling the coupling strength appropriately. 

In contrast, for high temperatures and for an Ohmic 
spectral density J(aj) = 2aa;e~"/"<^ [2g|, the Fourier 
transform of the correlation function becomes essentially 
fiat for large cutoff frequencies ojc with a limiting scaling 
7 = k^T2a « 7o « 7±n. Here, the Markovian approx- 
imation becomes exact as the bath correlation function 
becomes a Dirac-(5 distribution. The high-temperature 
stationary state is then just the completely mixed state 
with <^CT^/^/^) = and in essence, decoherence occurs 
more rapidly at higher temperatures. In this limit also 
further terms corresponding to the Lamb-shift are neg- 
ligible. For simplicity and to stress the drastic merit 
of feedback control, we will therefore in the following 
focus on this worst-case scenario for decoherence (i.e., 
the high-temperature and wide-band limit), where the 
Fourier transform of the bath correlation function is char- 
acterized by the single parameter 7. 



B. Quantum Zeno Limit 

Performing repeated measurements but neither allow- 
ing for feedback nor dissipation (A± = 0), the effective 



propagator ([7]) becomes 



Pzono(Ai) = - 



(1 1 








1 1 


\ 




g-iOAi 
g+iOAt 




g-is:2At 

g+iOAt 



(11) 



which simply leads to the Quantum Zcno effect [29|, |3C 
Repeated application of the above propagator yields 



^Zc„o(At) 



1 

cos""i(f7At)l 



7'zc„o(At) , (12) 



which when t = rtAi is kept constant while At — >■ and 
n — >■ oo reduces further to V^^^^{/S.t) — > 7'zcno(0). Obvi- 
ously, this operator preserves the pure eigenstates of the 
measurement operators \^) = -y= [|0) ± |1)], such that 
the system when initialized in these states will remain 
frozen. 



C. Zeno and dissipation 

Adding dissipation but without feedback (A± = A > 
and n± = n), the Quantum Zeno effect does not sta- 
bilize the coherences: According to Eq. ([T]), these decay 
between the measurements as exemplified by the decay- 
ing expectation values 



{<^t+At) - cos(r!At)e'''^*^'[3+™^(2'')l/2 (af) 
i'^t+At) = sm{nAt)e 



-7AtA^[3+cos(20)]/2 I x\ 



, (13) 



whilst (cr^_(_At) = 0- Note that the chosen measurement 
of (T^ at time t projects ( a^'^ \ to zero, which eventu- 
ally implies that only the expectation value of (erf) has 
a differential equation limit as At — >■ 0. 



sensitive to the ratio of the dampening constants A±: 
The stationary state is purified when A-j- 3> A_ (where 



— 7^ — 1) or A+ <C A- (where 



-1 



The 



weak dependence on the dissipation direction demon- 
strates that small control errors (e.g. caused by further 
decoherence channels) have only small effects on purifica- 
tion efhciency. When the dampening rates differ strongly 
(strong feedback) , the eigenstate of the measurement op- 
erator that has the smaller dampening rate becomes pu- 
rified. The applicability of the effective description of 
measurement and feedback control by a fixed-point iter- 
ation can also be checked by comparing with averaging 
over many different trajectories, see Fig. [TJ 




D. Feedback 



0,4 0,6 

time y t 

FIG. 1: (Color Online) Comparison of the effective evolution 
of the (erf) expectation value under feedback control for finite 
measurement intervals (symbols) with an average of 10^ , 10'' , 
and 10^ trajectories (thin dotted, dashed, and solid curves 
in lighter colors, respectively; single trajectories would decay 
at different pace from ±1 towards between measurements). 
In the continuous measurement limit (top bold curve), the 
final purity is maximal and would vanish without feedback. 
Parameters: A+ = 1.0, A_ = 5.0, fl/j = 5.0, e± = tt/2. 



Finally, with both dissipation and feedback, different 
conditioned Liouvillians are applied between the mea- 
surements. Using the Bloch sphere representation p ~ 
^[1 + r ■ cr] with |r| < 1 we can rephrase the fixed-point 
iteration for the density matrix as a fixed-point iteration 
for the expectation values of the Pauli matrices, see ap- 
pendix [Cj Of these, only one has a continuum limit when 
Ai^O 



(^D 



7 



{A^ [3 + cos(20_)] + A^ [3 -I- cos(26l+)]} (erf) 



+ 1 {xl [3 + cos(26l_)] - A^ [3 + cos(26l+)]} ,(14) 



whilst the other expectation values are continuously pro- 
jected to zero. First we see that under feedback, the sta- 
tionary solution (fj^) does not vanish despite the high- 
temperature limit. Remarkably, it is only weakly de- 
pendent on the interaction angles 9± but is much more 



Unless either A_ or A+ vanishes (which effectively in- 
troduces a decoherence- free subspace), the continuous 
measurement limit At — >■ is one necessary ingredient 
for purification: Then however, occurrence of purifica- 
tion is not specific to the chosen measurement direction. 

A similiar calculation can be performed with 
the Hamiltonian 1^ for a general measurement di- 
rection M± = i [1 ± n • cr] defined by n = 
(sin(6') cos(0), sin(6') sin(0), cos{9)), see appendixiDl Note 
that in this light, the Maxwell-demon type classical feed- 
back in Ref. |3l| appears as a special case where system 
Hamiltonian of a single- level quantum dot Hg = ed^d 
with fermionic operators and measurement operators for 
empty Me = 1 — d^d and filled Mp — d'^d dot states com- 
mute. For a general measurement direction one finds for 
our model that the purification direction coincides with 
it. The modulus of the Bloch vector of the stationary 



state does not depend on (/> and 



values 



E 



[xif{9,e+) + xu{e,0^)y 



(15) 



and only weakly on the remaining angles f{9, 9i) = 5 — 
cos(20i) — cos(20) [1 + 3cos(20i)] (the special case 9 — 
9i = corresponds to trivial purification in the direction 
of Hs). It is immediately evident that it approaches one 
(pure state) when one of the dampening coefficients A± 
vanishes or is much larger than the other. 



(^ 



xx\ 



(E^^) =+(S^^) 



X2 



X2 



'^R + 3Ab 



(18) 



see also appendix |F1 Similarly, the pure Bell state is 
fully characterized by the only non-vanishing expecta- 
tion values (E^^) = - (E'^^) = (S^^) = +1. Eq. HH]) 
corresponds to a stationary concurrence [33| and purity 
P = Tr{p2} of 



III. DISSIPATIVE ENTANGLING 

For two qubits, one would be interested in a measure- 
ment scheme that leads to a purified entangled state. 
A projective (non-demolition) measurement of the Bell 
state I^Bcii) = [|00) -I- |11)] /\/2 corresponds to a two- 
local measurement Bell (B) projection operator 



Mb 



I'J'f 



(* 



Boll 



4 ^ 



CTi cr 



1"2 



^1^2] -(16) 



Assuming the measurement only to recognize the Bell 
state (two outcomes), the remaining (R) projector fol- 
lows from the completeness relation as Mr = 1 — M-q ■ In 
the energy eigenbasis defined by the local system Hamil- 
tonian 



UJi u;2 



(17) 



the 16 X 16 superoperator equivalents of MbpM^ and 
M^ipM^ will not have a Block structure between coher- 
ences and populations, see appendix [E] For simplicity, 
we choose to change the strength of the interaction only 



^SB — A 



B/R (o-f -I- erf ) (Ki Y.k [hkbk + hlh\\ and to keep 

the bath invariant as in Eq. ^. The resulting Liou- 
villians in Born-Markov-secular approximation decouple 
populations and coherences in the system energy eigen- 
basis and simplify strongly in the high-temperature limit 
7tj = 7 (where also the Lamb-shift vanishes), see ap- 
pendix |E1 According to Eq. ([T]) , the effective propagator 



becomes 7^cff(At) 



,£BAt 



Ml 



oCnAt 



A4r and defines 



a fixed-point iteration for the density matrix. Just as for 
a single qubit one may parameterize the density matrix 



by Pauh matrices [3^ p(t) = T,ai3e{o 
with cr°/2 = 



^}' 



r°^l3 



{t)afa: 



I1/2 and r"" — 1/4. The generators 



ct"iT2 of the group 5J7(4) are trace orthogonal 
Tr|l]"^E"'^'| = 4:Saa'Si3i3', which allows to obtain a 
fixed-point iteration for the 15 nontrivial expectation val- 
ues of generalized Pauli matrices ( S" \ — 4r"'^ (i) . The 

stationary state of the fixed-point iteration (also identifi- 
able as the normalized eigenvector of the effective prop- 
agator in Eq. ([7]) with eigenvalue one) has in the con- 
tinuum limit At -^ the only non-vanishing stationary 



C = ^^ 



V + 3Ab 



'^Be(A^-3A|), P= ..T + '?.2 (19) 



(A| 



3A^ 



with Q{x) denoting the Heaviside step function. Thus, 
when the undesired parts of the density matrix are 
damped much stronger than the entangled parts Ar, ^ 
Ab, both concurrence and purity approach one, which 
demonstrates that the method may in principle create 
pure maximally entangled states. In contrast, without 
feedback (Ar = Ar), the stationary concurrence vanishes 
and the purity becomes that of a completely mixed two- 
qubit state. For finite measurement intervals At, the 
fidelity of the state preparation reduces in comparison 
to Eq. P^ . see Fig. [31 A perfect stabilization of the 
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FIG. 2: (Color Online) Effective evolution of the non- 
vanishing ( E( ) expectation values under feedback control 

for finite measurement intervals (symbols) and in the con- 
tinuum limit (solid curves). Parameters: Ab = 1, Ar = 5, 
(jJA^t — ujBAt = 0. A non-vanishing stationary concurrence 
(see inset extending Eq. (I19|l to finite measurement intervals, 
with contour lines ranging from 0.1 to 0.9 in steps of 0.1) 
requires short measurement intervals. 

Bell state therefore requires both the continuum limit 
where measurement, processing, and control are per- 
formed much faster than decoherence, and strong control 
Ar/Ab -)■ 00. 



IV. SUMMARY 

For a piecewise-constant and instantaneous control, it 
is possible to explicitly include the measurement pro- 
cess into a quantum feedback scheme for the evolution 
of an open quantum system. However, unlike standard 
quantum feedback control, an effective description by a 
master equation is not possible anymore. Instead, the 
evolution may be described by a fixed-point iteration of 
the density matrix. As a very intuitive outcome, purifi- 
cation of eigenstates of the measurement operators can 
be achieved provided that unwanted parts of the density 
matrix are damped with a significantly larger rate than 
the desired ones. By construction, further decoherence 
channels simply enter the scheme as control errors, for 
which we have found only a mild effect on the purifica- 
tion efficiency. The scheme is suitable for high tempera- 
tures, but generally the time interval at which measure- 
ments are performed has to be much smaller than the 
temperature-dependent decoherence time corresponding 
to the conditioned Lindblad evolutions. Therefore, even 
when measurement-induced collapse and control are in- 
stantaneous, the purification efficiency will eventually 
be limited by the finite signal processing speed thereby 



for classical electronic processing requiring decoherence 
times in the order of milliseconds. 

To achieve significant purification without a 
decoherence-free subspace, it is important to note 
that the limitation to Markovian and weakly-coupled 
systems used in the derivation of the Lindblad dissipators 
can in the continuum measurement limit be overcome 
by deriving the trace- and positivity-preserving map 
via coarse-graining [2J, [2^, where the coarse-graining 
timescale is naturally set by the detector sampling 
frequency At. A further interesting avenue of research 
is how the quantum version of the Jarzinsky equal- 
ity [33, [33 is modified by feedback and non-pro jective 
measurements that do not commute with the system 
Hamiltonian. 
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Appendix A: Jump measurement limit 

Suppose we consider a two-level system under the continuous action of a vaccuum reservoir and a measurement 
device that detects whether during time intervals At a quantum jump has occurred in the system (e.g. decay to the 
ground state via the emission of a photon) . We assume that the decay of the two- level system from its excited state 
|1) to its ground state |0) is detected with unit efficiency and occurs with probability 7At. A measurement of such 
quantum jumps requires a POVM description, since these are non-projective. In addition, the measurement operators 
for the outcomes click (c) and no-click (nc) will depend on At 



Mc(At) = V7At|0)(l| , M„e(Ai) = v/l-7At|l)(l| + |0)(0| , (Al) 

which automatically obeys the completeness relation M^{At)Mc{At) + Ml^{At)AInc{A.t) — 1. They act on the density 
matrix as 

M,{At)pMl{At) = 7Aipn |0) (0| , 



M^,{At)pMl{At) = Poo |0) (0| + (1 - jAt)pn |1> (1| + ^1 - lAtpoi |0) (1| + \/l - -fAtpw |1) (0| . (A2) 
In the following, we will need their action for small At and therefore state 
MciO)p = Afe(0)pMt(0) = O, 

Mnc{0)p = Afnc(0)M4(0)=P, 

M'Mp - ^[Mc(At)pAft(Ai)]^^^p-7Pn|0)(0|, 
MU0)P = ^ [M,{At)pMl{At)]^^^^ = -7P11 11) (1| - Jpoi |0) (1| - |pio |1> (0| . (A3) 

The feedback scheme can now be defined by performing an instantaneous unitary transformation Uc (experimentally 
approximated by using a (5-kick like pulse on the system Hamiltonian) of the density matrix whenever a detector click 
has occured. In addition, the system Hamiltonian (and possible further unmonitored reservoirs) is included in the 
Lindblad superoperator Co- Then, the effective propagator becomes 

ViAt) = e^^^'e'^-MciAt) + e^°^'Mnc{At) , (A4) 

where e'^" p = UcpUl and correspondence to Eq. (jH]) is established by putting Bc{At) = e'^o'^'^e'^" and Bnc{At) — e'-'°^^ . 
For small At (continuum limit), this allows to define an effective master equation under jump detection and unitary 
control 

p = hm P^'^^')-P^^ ^ linr J- [V{At) - 1] p{t) 
^ At^a At At^o At^ ^ ' ^^ ' 

= lim -i- {[e'^'MciO) + MnciO) - 1] + At [Coe^'Mc{0) + CoMnc{0) + e^^M'M + M^O)]} p{t) 

At-)-0 At 

= [+£o + e^^M'M + M[,M\ pit) - Cf^p{t)^Ci^[p{t)] , (A5) 

where we have already used the superoperator equivalents of Eq. (jA3p . At the second line, a necessary ingredient for 
an effective master equation description becomes obvious: To remove the first term in square brackets, it is formally 
required that '}2,^Bm{^)M.m{^) = 1- Here, this is fulfilled since infinitely short measurements do not affect the 
density matrix. Finally, also inserting the derivatives of the measurement superoperators makes the Liouvillian in the 
effective master equation picture explicit 



Cih[p] = Co[p] +7 
= C,q[p] +7 



C/c|0)(l|p|l)(0|f/t-i{|l)(l|,p} 
L,pLl-]^{LlL,,p) 



(A6) 



with Lc — Uc |0) (1|, which obviously is a master equation of Lindblad type. 



Appendix B: Single-qubit Liouvillian 



The Born-, Markov-, and secular approximations lead for a non-degenerate system Hamiltonian generally to a 
decoupling of the evolution equations for populations and coherences in the system energy eigenbasis. For the single 
qubit with fi > 0, non-degeneracy is obviously fulfilled. Starting from Eq. (|9]) in the main text, the Liouvillians act 
on the density matrix elements pij = {i\p |j) as 

Poo = -A| sin^(6'±)7+j2Poo + A| sin^(6'±)7_opii , 

Pii = +A±sin^(6'±)7+nPoo - A±sin^(6'±)7_opii , 

Poi - [-i^ - A| siv?{6±) {a+n - ct_o) /2 - A| sixi^{6±) (7+^ - 7_o) /2 - 2A| cos2(0±)7o] poi , 

pio = [+iO + A^ siv?{6±) {a+n - c7_o) /2 - A^ sir?{6±) (7+n - 1-n) /2 - 2A^ cos'{e±)-io\ pw , 

where 71^ is the even (real-valued) Fourier transform of the bath correlation function 



C(t) = (e 



-iHbt 



Be 



-iHsT 



B) 



E 

kk' 



hkbke 



hlb[e 



t^+it- 



hk'hk' +hl,b 



CXJ 

J2 \hkf [e-'^'^ [1 + neK)] + e+^^-^neK)] = -^ / J(w) [ 

7, "^ 



(Bl) 



nBiuj)] + e+'"^^nBiuj)]duj 



+ CX3 



— / JH[l + nBH]e" 



'duj 



(B2) 



as defined in the main text and similarly the imaginary quantity 



(^ 



-iHbt 



Be 



-iHsT 



B) 



+ ILJT 



sgn (r) dr ~ —V 



-du) 



(B3) 



denotes its odd Fourier transform, which may also be obtained from the even one via a Cauchy principal value integral. 
The odd Fourier transform is associated with the Lamb-shift terms that account for an energy renormalization of 
the qubit. Note that in the high-temperature and wide-band limits that we focus on, 71^ is approximately flat, such 
that cr^j ~ 0. Then, the stationary state of both Liouvillians is the completely mixed one with poo = Pii = 1/2 
and poi = Pio = 0. However, tJtj is purely imaginary in any case, which for the evolution of the absolute square of 
coherences implies Eq. (|10p in the main text. Similarly, we obtain thermalization from the evolution of populations 
by invoking the KMS condition. 

Appendix C: Bloch sphere iteration 

The Bloch sphere representation of the density matrix p — ^[1 + r ■ cr] implies {a") ~ r" or alternatively for the 
density vector the representation 



P = ( ^ (1 + i^l) , ^ (1 - ('^^» , I i{<^l i (^^» , I ((^^) + i ('^^» ) , 



(CI) 



which can be inserted in the fixed-point iteration for the density matrix p{t + At) = Vcs{A.t)p{t) to yield iteration 
equations for the expectation values of Pauli matrices 

cos(r2At) 



/„x \ _ ^^^V^^J f -fAtXl[3+cos{2e+)]/2 ^-7AtA^ [3+cos(2e_)]/2 

\'^t+At/ ~ 2 \ 



I COS(iJAtj / _^AtA^[3+cos(2e+)]/2 ^g-7AtA?.[3-f^cos(2e_)]/2\ ^^x'^ 



i'^t+At) 



sm{nAt) 



-7AtA|^[3+cos(2e+)]/2 _ -7AtA^ [3+cos(2e_ )]/2 



I SmjiZAtj / -yAtA^f3+cos(2e^)]/2 ^ ^-7AtA!. f3+cos(29_)]/2\ ^^x-^ ^ 



i^t+At) = 0. 



(C2) 



The ((T^)-expectation values does not recover from zero after the first projection since the chosen system Hamiltonian 
was proportional to a^ , which only generates rotations between a^ and a^ . These equations recover the no- feedback 
limit (A-t = A and B± = ff) disscussed and the continuum limit in Eq. (|14p in the main text. 



Appendix D: Most general measurement direction 

For more general projective measurements M± = i [1 ± n • <t] = Mj., their action on the density vector in ordering 
P= (poo,/3ii,Poi,/5io)'^ is given by 



1 



M^ = 4 



M- 



1 

4 



/ 4cos4(f) sin2(0) +4e+'^sin(f)cos3(|) +4e-''^sin(f)cos3(|) \ 

sin2(6l) 4sin''(f) +4e+''^sin3 (f ) cos (|) +4e-''^sin3 (|) cos (|) 

+4e-'^sin(|)cos3(|) +4e-''^sin3 (|) cos (|) fii^{e) e-^'^sii?{d) 

V+4e+'^sin(|)cos3(|) +4e+''^sin=^ (f ) cos (f ) e+^'* sin\e) sm''{9) J 

( 4sin4 (f) sin2(0) ~4e+'*sin3 (f) cos (|) -4e-'*sin3 (f) cos (|) \ 

sin2(0) 4cos4(f) -4e+'^sin(|)cos3(|) -4e-''^sin (|) cos^ (|) ,. 

-4e-''^sin3(|)cos(|) -4e-'"*sin (|) cos^ (|) sin2(6') e-'^'"f' srn^ {9) 

V-4e+'^sin3(|)cos(|) -4e+''^sin (|) cos^ (f ) e+2'^sin2(0) siu^iO) j 



;di) 



One can easily check the orthogonal projector properties M'^ = M.± and M.±M^ = 0. In addition, the special case 
of Eq. (4) in the main text is obtained by putting 6 — 7r/2 and 4> = 0. Inserting the corresponding effective propagator 
into the fixed-point iteration for the Pauli matrix expectation values one now observes that all expectation values 
couple to each other. Solving for the stationary state yields Eq. P3|) of the main text. Naturally, Eq. ((T5|) is also 
compatible with the stationary state of Eq. dTi|) when — -k jl and = 0. 



Appendix E: Two-qubit Liouvillian 

In the energy eigenbasis of the two qubits |00), |01), |10), |11), the Born-, Markov-, and secular approximations again 
lead to a decoupled evolution of populations and coherences. The populations {pij,ki = {ij\ P 1^^)) evolve according to 

Poo, 00 = ^('^b/rT+cji + •^b/r7+<.j2)poo,oo + AByf{7_cj2Poi,oi + '^b/r,7-wiPio,io , 

^01,01 = +AB/r{7+cj2Poo,oo - ('^b/rT+wi + Ab/r7-w2)poi,oi + ^b/b^I-^iP^i^^'^ > 

Pio,io = +Ab/r, 7+^1^00,00 - (^b/r7-wi + Ab/r7+cu2)pio,io + '^b/r7-!.j2Pii,ii ) 

Pii,ii = +Ab/r,7+wi/3oi,oi + Ab/r7+w2Pio,io - ('^b/r7-wi + Ab/r7-w2)pii,ii (El) 



and the coherences decay independently (omitted for brevity). By invoking the KMS condition ^^^ = e^'^"7+tj it 
is immediately obvious that the stationary state is the thermalized one. In the high-temperature limit considered in 
the paper, the stationary state of either stationary Liouvillian is just the completely mixed state. In this limit, the 



coherences decay according to 



10 



Poo, 01 = 

Poo, 10 = 

Poo, 11 = 

Poi,oo = 

Poi,io ~ 

P01,ll = 

PlO,00 = 

PlO,01 = 

PlO,ll = 

Pll,00 = 

Pll,01 = 

Pll,10 = 



-\UJ2 



- 2A|/j,7^ 



Poo, 01 + -^b/r7Pio,ii , 



-icji - 2Ab/p7J Poo, 10 + '^b/rTPoi.h , 
' 2A|/r7J poo,ii , 

\2 



-icji 



\LlI\ — iW2 

-iu;2 - 2Ab/j^7J poi,oo + Ab/j^7Pii,io , 
+ iw2 - 2A|/j^7J /9oi,io , 
-icji - 2A|/j^7J poi,ii + Ab/r7Poo,io , 
-icji - 2Ab/j^7J Pio,oo + Ab/j^7Pii,oi , 
iw2 - 2A|/j^7J pio,oi 



\-\UJ\ 



01 , 



— \UJ2 — 2A: 



3/r7J Pio,ii + Ab/r7Poo,oi , 
hicji + iw2 - 2Ag/p7J pii^oo , 
hicji - 2A^/p(7J pii,oi + A|/r7Pio,oo , 
- 2Ab/r7J Pii,io + Ab/r7Poi,oo , 



\-\UJ2 



(E2) 



where coherences may couple among themselves provided they belong to superpositions of states with equal energy 
differences. Inserting 7^,^ — >■ 7 also in the equations for the populations yields together with the sparse 16 x 16 
measurement superoperators 
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(E4) 



acting from the left on density vectors arranged as 

( Poo,ooj Poifii, Pwao, Pii, Hi Poo, 01, Poo, 10: Poo,iij Poi,ooj Poi,iOj Poi,iij Pio,ooj Pio,oij Pio,ii, Pii.oo, Pii.oi, Pii.io ) , the ef- 
fective propagator for the fixed-point iteration of the density matrix. 



Appendix F: Stationary state of two-qubit fixed-point iteration 



Mapping the fixed-point iteration for the density matrix to the expectation values of (S"^), one obtains lengthy 
expressions for the 15 nontrivial expectation values: (cf), (o-|), (cr|), {(jf), (ti), (cf), (erf erf), {(^1(^2)^ i'^i'^i): 
(cTj'crf), {cr\(72), (erf erf), {(^{(72), {afcr"^), and (erf erf). These can be solved for the stationary state. Alternatively, 



the stationary state of the iteration can be obtained directly by looking at the trace-normalized eigenvector of the 
effective propagator with eigenvalue one. For finite At, it is characterized by vanishing expectation values of all local 
operators. The non-vanishing stationary expectation values read 



'l'^2/ 



'l'^2/ 



(^I'^s; 



lo-iCr2; 



4g3(AB-HAR) sinh [Ab - Ar] sinh [2Ar] sin(f7) 



2Ae 



e2AR)cos(f^) 



gOAR -(_ g4AB+2AR (3 _ 4e4AR) _|_ (e2AB _|_ gSAR^ (-2e2AB+4AR _ g 
4g3(AB+AR) gjjjjj [^g _ ^j^J gjj^j^ [2Ar] COs(r2) 
eOAR _^ g4AB-(-2AR (3 _ 4g4AR) _^ (g2AB + e^An) (2e2AB+4AR _ g2AB _ gZAR) cos(f2) 
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. g4AB-|-2AR ^4g 



4A„ 



3) 



o4Ab 



o4An 



2e4AB+4AR _ 2e2AB+2AR (g. 
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1)] cos(r2) ' 



'1"2/ 



(Fl) 



where we have used the dimensionless quantities ft = (wi + W2)At, Ab = 7AiA|, and Ar = 7A<Ap. In the limit 
At — >■ 0, this recovers Eq. (11) in the main text. Even for finite At however, the stationary density matrix is just 



given by three independent parameters ai 



[at a. 



1"2/ 



0:2 



rffff) 



(a^a^), and a^ = (erfcr|), for which 



one may also express concurrence and purity analytically which generalizes Eq. (1191) in the main text (not shown). 



